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We propose a scheme for parametric amplification and phase conjugation of an atomic Bose- 
Einstein condensate (BEC) via stimulated dissociation of a BEC of molecular dimers consisting 
of bosonic atoms. This can potentially be realized via coherent Raman transitions or using a 
magnetic Feshbach resonance. We show that the interaction of a small incoming atomic BEC with 
a (stationary) molecular BEC can produce two counterpropagating atomic beams - an amplified 
atomic BEC and its phase-conjugate or "time-reversed" replica. The two beams can possess strong 
quantum correlation in the relative particle number, with squeezed number-difference fiuctuations. 
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PACS numbers: 03.75.Kk, 03.65.Yz, 42.50.-p 

The fascinating experimental progress in controlling 
ultracold atomic gases has now resulted in the creation of 
quantum-degenerate samples of ultracold molecules and 
molecular Bose-Einstein condensates (BEC). The first 
step toward seeing molecular condensation was under- 
taken in transient experiments with a BEC of ^^Rb atoms 
[l| , in which interference measurements were indicative of 
a small molecular condensate formation. More recent ex- 
periments with ^^•^Cs, ^^Rb, and ^■^Na ", as well as with 
degenerate Fermi gases of ^°K and ®Li atoms |3, have 
produced even larger samples of ultracold molecules, in- 
cluding molecular condensates. 

All these experiments have employed conversion of 
atom pairs into weakly bound molecular dimers in the 
vicinity of a magnetically tunable Feshbach resonance. 
This technique appears to be more successful at present 
than Raman photoassociation |j, |5(. Unambiguous 
claims on the observation of equilibrium molecular con- 
densates have only been made for the case of molecules 
consisting of fermionic atoms. These molecules have 
much longer lifetimes because of Pauli blocking, which 
suppresses their decay due to inelastic collisions |g. De- 
spite this, the production of Bose condensed dimers com- 
posed of bosonic atoms should not be out of reach using 
current experimental techniques. 

The purely bosonic case is particularly interesting since 
the underlying dynamics in the atom-molecule conversion 
can take advantage of coherence and bosonic stimulation 
with respect to both the atomic and molecular species. 
In practical terms, this could lead to wider range of tech- 
niques for coherent quantum control of ultracold quan- 
tum gases, including the possibility of superchemistry [5J, 
matter-wave amplification with atom-molecular "laser" 
beams, as well as various implementations of nonlinear 
atom optics and quantum atoms optics. 

Here, we propose a scheme for nonlinear, coherent 
matter-wave amplification and phase conjugation, which 
at the same time is capable of producing pair-correlated 
or number squeezed atomic beams. The scheme relies 



on the process of stimulated dissociation of a BEC of 
molecular dimers, taking place in the presence of an in- 
jected signal - a small incoming atomic condensate. The 
resulting output consists of an amplified signal and its 
"time-reversed" or phase-conjugate replica [71 . 

This is analogous to the nonlinear optical process of 
parametric amplification of a "seed" signal, using fre- 
quency conversion of photons in a quadratically nonlin- 
ear media [S|. In the matter- wave parametric amplifica- 
tion, the coupling that takes the role of the "quadratic 
nonlinearity" and converts molecules into atom pairs and 
vice versa can be realized via coherent Raman transitions 
or using magnetic Feshbach resonances. Either of these 
mechanisms have certain advantages and disadvantages 
from the practical point of view, however, the essential 
physics can be modeled via an effective quantum-field 
theory, which is identical in the two cases |3, U^ llj • 

The advantage of the present scheme compared to 
spontaneous molecule dissociation [1^ [ij, [ij, Ha, llq| . 
with no atomic condensate present initially, is that the 
desired quantum effects with mesoscopic atom numbers 
can be achieved on much shorter time scales. As a result, 
the disruptive effect of inelastic collisions on molecule life- 
times can be bypassed, thus making the present scheme 
much more robust for practical implementation. In ad- 
dition, the scheme has an advantage compared to earlier 
related proposals |l7| using four-wave mixing |l8| in that 
it is less susceptible to phase noise in the "pump" BEC 
- once its depletion is taken into account. 

We start by considering the Heisenberg equations of 
motion for the coupled atomic-molecular system in a one- 
dimensional (ID) environment, which in a rotating frame 
and in a dimensionless form are given by [9|, ll3J 
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Here, i^ii^^r) = v^*i(s,i) and ip2{^,T) = %/d^*2(a;,i) 
represent the dimensionless atomic and molecular field 

operators (such that / ^|^i ) gives the ID linear density), 

with bosonic commutation relations ['4'i{£,,T),iljJ£^' ,t)] — 
SijS{^ — £_'), where S, — x/do is the dimensionless coordi- 
nate, T = t/to is the dimensionless time, do is a length 
scale, and to — 2mid^lh is a time scale. In addition, 
TOi and 7Ti2 — 2mi are the atomic and molecular masses, 
6 = Ato is the dimensionless detuning (with 2hA giving 
the overall energy mismatch between a pair of free atoms 
and a molecule), k — xto/Vdo is the dimensionless atom- 
molecule coupling (with x being the respective ID cou- 
pling), and U22 = U22to/do is the dimensionless coupling 
for molecule-molecule s-wave scatterin g int eraction (with 
U22 being the respective ID coupling [l^ll^). Finally, 
^2(0 — — "22«2(1 ~ C^/^o) describes the axial trapping 
potential for the molecules which we assume is harmonic, 

where n^ — (il>l{0, O)'02(O, 0) ) is the dimensionless initial 
peak density of the molecular condensate and ^0 is the 
respective Thomas-Fermi radius. 

We assume a highly elongated cigar-shaped trap ge- 
ometry such that the system can be modelled by a ID 
theory []j|, |20|, as above. The atoms are assumed to 
be untrapped axially, while confined transversely. The 
molecules are trapped both transversely and axially. 

In the case of Raman photodissociation, the energy 
mismatch is given by 2HA = 2Ei — E2 — Sw, where w is 
the frequency difference between the two Raman lasers, 
while 2Ei and E2 refer, respectively, to the energy of 
free-atom pairs in the dissociation limit and the energy 
of the bound molecular state. 

In the case of a Feshbach resonance [l5|, [ig, 2/iA 
gives the energy mismatch achieved upon switching on 
the atom- molecule coupling, i.e., upon a rapid cross- 
ing of the magnetic field through the resonance, which 
"brings" the (initially stable) molecular level above the 
atomic dissociation limit. This corresponds to having 
A < 0, and results in dissociation of molecules into atom 
pairs such that the potential energy 2?i|A| is converted 
into the atomic kinetic energy [2?i|A| ~ 2?i^fc^/(2mi)] for 
selected phase-matched modes with opposite momenta 
around ±fco = yj2mi\A\/h. This strategy of a fast ramp 
through the resonance is very similar to the one real- 
ized experimentally in Ref. 0, where a quasi-mono- 
energetic spherical wave of atoms was created. A real- 
ization of this experiment in a one-dimensional environ- 
ment would correspond to the conditions of the current 
proposal. 

In our model, the interaction terms due to atom-atom 
and atom-molecule s-wave scattering processes are ne- 
glected on the grounds that we only consider large ab- 
solute values of the detuning |A| and restrict ourselves 
to short dissociation times. As a result, the number of 
atoms produced during dissociation remains small, and 
the mean-field or phase diffusion terms due to atom- 



atom and atom-molecule scatterings remain negligible 
compared to |A|. 

Before analyzing the quantum dynamics of the sys- 
tem described by Eqs. 1^, we first consider a simplified 
model that has an analytic solution. This corresponds 
to an undepleted, uniform molecular condensate in a co- 
herent state at density n2, in which case the molecu- 
lar field amplitude (which we assume is real) can be ab- 



sorbed into an effective coupling g = ny 



Expanding 



-01 (^,t) in terms of single-mode annihilation operators, 

V-'i(Cj''') = X]g ^g(''')e*''^/v^, where q — kdo is the di- 
mensionless momentum, I is the quantization length, and 
the operators dq satisfy the usual commutation relations 
[dqir), a , (r)] = Sq^q' , we obtain the following Heisenberg 
equations of motion: 



-i [q^ + d]dq + gdLq 



ddq/dr 

ddLq/dr — i[q'^ + S] aL^ -t- gdq. (2) 

These have the following well-known solution [ij, 
113: dqir) = Aq{T)aq{0) + Bq{T)dlq{0) and algir) = 
BqiT)dq{0) + A*{T)dl^iO), whcrc 



Aq{T) = cosh (gqT) - iXq Sllih (gqT) /gq, 
Bq{r) = g sinh {gqT) /gq, 



(3) 



and \Aq\ — B^ = 1. Here, the parameter Xq = q-' + 6 
can be identified with an effective phase mismatch, while 
gq = [g^ — A^)^/^ is the gain coefficient, which - if real - 
causes a growing output in mode (7, while ~ if imaginary 
- leads to oscillations. 

To study parametric amplification in the atomic field, 
we consider initial conditions where all momentum com- 
ponents with q > Q are initially in a coherent state 
0,^(0) \aq) — Uq \aq)^ whcrc Uq are the corresponding am- 
plitudes, while all negative momentum components are 
initially in the vacuum state d_g(0) |0) = 0. 

Next, we introduce particle number operators N^ — 
5I)(7>o ^\^q ^^d N- — X]g>o OL-qO'-q corrcspondiug to the 
total number of atoms with positive and negative mo- 
menta, respectively, and find that the respective average 
numbers are given by: 



^^+i-))=Y.q,,[Bl{r){l 



\aa 



{N-{t)) 
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(4) 
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For (5 < 0, the function B^ has two distinct global max- 
ima at ±(7o = vPl corresponding to a zero effective 
phase mismatch, A^ — 0. In the expression for {N^{t)), 
the three terms under the sum are identified as the am- 
plified contribution of the vacuum noise in the mode 
q, B'^{t), the amplified coherent component of the in- 
put, i?^(T)|agp, and the coherent input component it- 
self, la^P, while {N-{t)) consists of the amplified vac- 
uum noise and the phase conjugate of the input. 



To quantify the correlation and relative number 
squeezing between 7V_|_ and -/V_, we consider the nor- 
malized variance V{t) of the particle number difference 

N+{t) — iV_ (r). In normally ordered form, V{t) is given 
by 

Vir) = 1 + (: [A(iV+ - iV_)]2 :)/((7V+) + (7V_)). (6) 

Here, AX = X ~ {X), and V{t) < 1 implies squeezing 
of fluctuations below the coherent level, which is due to 
strong quantum correlation between the particle numbers 
in 7V+(t) and N_{t). 

Calculating the quantities (: (7V+^_)^ :) and (N^N^) 
in Eq. © and assuming that jckgp :§> 1 gives the following 
approximate result for the variance: 



V{t) o, 1 - 



E,>o[2i?,2M + l]KI 



(7) 



As we see, the degree of squeezing depends on the mag- 
nitude of the amplification factor B'^{t), and for strong 
amplification, B'^{t) 3> 1, one can obtain almost perfect 
(100%) squeezing, V{t) ~ 0. 

We now turn to the exact quantum dynamical simula- 
tion of the nonuniform system, Eqs. J^). Here, we take 
into account molecular field depletion, molecule-molecule 
s-wave scattering, and we include possible (linear) losses 
of atoms and molecules, occurring at rates 71 and 72, re- 
spectively. The simulation is done via numerical solution 
of the stochastic (c-number) differential equations 13] in 
the positive-P representation [22| 
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- (71 + iS)^pl + K'ijj2'ip^ + VK^hm, 

2 - [72 + «V'2(0 + «U22'0^'02] -02 



"tiPI + V-iU22'ip2m, 



(8) 



together with the equations for the "conjugate" fields 
ipi2, having noise terms ri^2- Apart from the new loss 
terms, these equations are equivalent to Eqs. (0, where 
ipi and ip^ are independent complex stochastic fields cor- 
responding, respectively, to the operators ij^i and tpj, 
while rji and rj^ {i = 1,2) are four real independent 
(5-correlated Gaussian noises with {r]i{^,T)rij{£_' ,t')) — 

Figure ^ represents an example of a quantum dynami- 
cal simulation of Eqs. ((HJ illustrating parametric amplifi- 
cation and phase conjugation of an incident atomic BEC 
with a center-of-mass momentum qg = ^/\6\■ Here, the 
top frame corresponds to the initial condition of a stable 
molecular condensate (shown in the middle) in a coher- 
ent state with a Thomas-Fermi density profile. The small 
incoming atomic BEC (shown on the left and moving to 
the right) is also assumed to be in a coherent state ini- 
tially, and since we neglect the atom-atom interactions 




position 

FIG. 1: Atomic and molecular (shown in the middle) density 
profiles, illustrating parametric amplification and phase con- 
jugation of an incoming atomic BEC (top frame, left) via 
stimulated molecule dissociation [23 • The results are ob- 
tained using 5,000 stochastic trajectory averages for k — 84, 
U22 = 1.8, S = -4.9 X 10*, 71 = 0.5, and 72 = 0. The disso- 
ciation coupling is switched on at ri = 0.01 for a duration of 
At = 4 X 10~*, while the total time window is r/ = 0.0204. 



in this low-density regime, we assume a Gaussian den- 
sity profile. The dissociation coupling k is invoked at 
time Ti when the atomic cloud is aligned with the molec- 
ular BEC (frame 5 from top). It is kept switched on for 
a short duration At such that the amplified and the "re- 
flected" output beams have densities comparable to that 
of the input beam and can be seen on the same graph. 

The next set of simulations is carried out for a more 
realistic set of parameter values than in Fig.^ We use a 
longer dissociation time and a larger coupling k - to re- 
sult in a large amplification factor and hence a strongly 
correlated output with squeezing in the particle num- 
ber difference. For simplicity, the simulation starts from 
T = Ti = when the incoming atomic cloud is already 
aligned with the molecular condensate and we invoke 
the atom-molecule coupling k. The dissociation is then 
stopped (k — 0) at T — T2, and we continue the dynami- 
cal evolution of the atomic field in free space (in ID) to 
allow spatial separation of the modes with positive and 
negative momenta during the time interval from T2 to tj. 

For spatially separated components, we can introduce 
a pair of particle number operators 



1/2(0) 



',U 



(9) 



N+i-){r)= / ^li^,r)MC,r)dC 

Jo{-l/2) 

Next, we define the normalized variance V{t) of the par- 
ticle number difference as in Eq. ^ and evaluate nu- 
merically the relevant averages, using the standard cor- 
respondence between the normally-ordered operator mo- 
ments and the c-number stochastic averages |22l |. 

Figure [3 shows the results for the atomic density dis- 
tribution ni{£_,Tf) — {'4'i{^,Tf)ipi{^,Tf)) at final time 
T — Tf and the variance V{t) as a function of t. In 
this simulation, the initial total number of molecules is 
415, the initial number of atoms in the incoming BEC is 




FIG. 2: (a) Final atomic density profile »^l(5,T■/) and (b) the 
variance V{t) as a function of time r. Here, the simulation 
(with 40,000 stochastic trajectory averages) starts at ri, the 
duration of dissociation is At = 8 x 10"'*, and Tf — 1.9 x 
10~^. Other parameter values are k — 297, U22 = 0.068, 
6 = -5.52 X 10'', 71 = 72 = 3, ng = 89.3, and Co = 3.49 H^. 



21, vifhereas the final number of atoms in the two output 
beams is ~ 60 and 78. The squeezing in the particle num- 
ber difference at t = tj is about 75% [V{Tf) ~ 0.25], and 
it is achieved on a time scale of At = T2 — ti = 8 x ID""* 
corresponding to 12.7 ms, with the parameter values used 
here |24J . This is much shorter than the time scale re- 
quired in spontaneous dissociation 13] to reach the same 
total number of atoms. At the same time, the degree of 
squeezing is still rather high. 

Thus, the disruptive effect of molecule losses due to 
inelastic collisions can be reduced in the present scheme. 
The reason for the shorter time scales required here is 
that the process of dissociation in the presence of a "seed" 
atomic BEC begins in the stimulated regime, with ex- 
ponentially growing output. In the case of spontaneous 
dissociation, on the other hand, the initial dynamics is 
in the spontaneous regime and the system spends a rela- 
tively long time here before the bosonic stimulation into 
the phase-matched modes becomes dominant. 

To summarize, we have analyzed the process of stim- 
ulated dissociation of a condensate of molecular dimers 
in the presence of a small incoming atomic BEC. This 
results in parametric amplification of the input BEC to- 
gether with generation of its phase-conjugate replica, 
propagating in the opposite direction. The two out- 
put beams are strongly correlated in the particle num- 
ber and have squeezed number-difference fluctuations. 
The squeezing with a mesoscopic total number of atoms 
can be achieved on much shorter time scales than in the 
case of spontaneous dissociation. This makes the present 
scheme more feasible for practical implementation, us- 
ing short-lived molecular condensates. In addition, the 
scheme provides a range of opportunities for coherent 
quantum control of ultracold quantum gases, including 
applications of nonlinear and quantum atom optics. 
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